source("F:/msaR.R")
Z<-read.table("clipboard",header = T)#????????
Z1=msa.X(Z)
View(Z1)
Y=scale(Z1,center=F,scale=T) 
View(Y)
D<-dist(Y,diag = T);D
plot(hclust(D,method = "single"))
plot(hclust(D,method = "ward.D2"))
H=hclust(D,method = "ward.D2")
plot(H);rect.hclust(H,4)
P=hclust(D,method = "median")
plot(P);rect.hclust(H,4)
C=hclust(D,method = "single")
plot(C);rect.hclust(H,4)
u=hclust(D,method = "average")
plot(u);rect.hclust(H,4)
n=hclust(D,method = "centroid")
plot(n);rect.hclust(H,4)